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O '■ ABSTRACT 

(~| . We present first INTEGRAL and XMM-Newton observations of a Seyfert galaxy, the type 

Ph' 2 AGN NGC 4388. Several INTEGRAL observations performed in 2003 allow us to study the 

spectrum in the 20 - 300 keV range. In addition two XMM-Newton observations give detailed 
^ ' insight into the 0.2 - 10 keV emission. The measurements presented here and comparison with 

C/2 , previous observations by BeppoSAX, SIGMA and CGRO/OSSFi show that the overall spectrum 

J^ ' from soft X-rays up to the gamma-rays can be described by a highly absorbed {Nh ~ 2.7 x 

10^^ cm~^) and variable non-thermal component in addition to constant non-absorbed thermal 
emission (T ~ 0.8 keV) of low abundance {Z ^ 7%Zq), plus a constant Fe Ka and K/3 line. 
, ; , The hard X-ray component is well described by a simple power law with a mean photon index 

C^ ■ of r = 1.7. During the INTEGRAL observations the 20 - 100 keV flux increased by a factor of 

1.4. The analysis of XMM-Newton data implies that the emission below 3 keV is decoupled from 
the AGN and probably due to extended emission as seen in Chandra observations. The constant 
iron line emission is apparently also decoupled from the direct emission of the central engine and 
likely to be generated in the obscuring material, e.g. in the molecular torus. 

Subject headings: galaxies: active — galaxies: individual (NGC 4388) — gamma rays: observations — 
X-rays: galaxies — galaxies: Seyfert 

1. Introduction 
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the broad line region, causing the characteristic 
Seyfert 2 optical spectra with narrow emission 
lines, but also absorbs most of the low energy X- 
ray emission. As the soft X-ray region is very sen- 
sitive to absorption in the line of sight, the study of 
the hard X-ray emission above ~ 3 keV is a power- 
ful tool to investigate the central engine of Seyfert 
2 (Ghisellini, Haardt & Matt 1994). In addition 
to this the variation of the Fe Ka line seems to 
vary differently with respect to the continuum in 
several cases (Iwasawa et al. 2003; Fabian 1977; 
Malzac & Petrucci 2002), which might be caused 
by the interaction of the core emission with the 
absorbing and reflecting material. 

At higher energies absorption becomes less 
important revealing the emission from the close 
environment of the central black hole. This 
energy range could be dominated by comp- 
tonized hard X-ray radiation from a hot (~ 
100 keV) plasma presumably forming an ac- 
cretion disk corona (Haardt et al. 1994), and 
Compton reflection of these hard X-rays on 
a cool accretion disk (George & Fabian 1991; 
Magdziarz & Zdziarski 1995). In Seyfert 2 galax- 
ies, reflection of the central radiation on the molec- 
ular torus often leads to an additional, possi- 
bly dominant, reflection component (Ghisellini, 
Haardt & Matt 1994). 

NGC 4388 (z = 0.00842; Philhps & Ma- 
lin 1982) is a good probe to test these sce- 
narios, as it is one of the brightest Seyfert 2 
galaxies at hard X-rays. The energy range up 
to ~ lOkeV has been studied with X-ray mis- 
sions hke ROSAT (AntoneUi et al. 1997), ASCA 
(Iwasawa et al. 1997; Forster et al. 1999), Bep- 
poS AX (Risaliti 2002), i?XT£; (Sazonov & Revnivt- 
sev 2004), and lately with Chandra (Iwasawa et 
al. 2003). Spectral studies of NGC 4388 of 
the hardest X-ray energies up to several hun- 
dreds of kcV have been performed by SIGMA 
(Paul et al. 1991), by the OSSE experiment (John- 
son et al. 1993) on the Compton Caninia Ray 
Observatory (CGRO; Gehrels et al. 1993), and up 
to ~ 150 keV by the Phoswich Detector System 
(PDS; Frontera et al. 1997) on-board BeppoSAX 
(Butler & Scarsi 1990; Boella et al. 1997). The 
two BeppoSAX observations in 1999 and 2000 
showed a high energy spectrum with a power law 
of r = 1.6 and 1.5, respectively. Due to the energy 
range covered by the PDS it was only possible to 



give a lower limit for a possible high energy cut- 
off {Ec > 109 keV). The INTEGRAL mission 
(Winkler et al. 2003), launched in October 2002, 
offers the unique opportunity to study the entire 
spectrum from 3 keV up to several MeV simulta- 
neously. NGC 4388 was the first Seyfert 2 galaxy 
detected by INTEGRAL in January 2003 (Beck- 
mann et al. 2004). In this paper we present data 
from several INTEGRAL observations and com- 
pare the results with previous SIGMA, OSSE, 
and BeppoSAX measurements. With a detec- 
tion significance of 10.9a and 12.6cr in the IN- 
TEGRAL/lSGRl energy bands 15 - 40kcV and 
40 — 100 keV, respectively, the data also allow us 
to study spectral variations between the 2003 ob- 
servations. In addition XMM-Newton data are 
presented to study the connection with soft X- 
rays. 

2. Simultaneous X-ray and gamma-ray ob- 
servations 

NGC 4388 was detected by INTEGRAL in the 
course of the 3C273 observation (Courvoisier et 
al. 2003a). Though NGC 4388 hes 10.4 degrees 
north of the quasar 3C273, valuable data for the 
Seyfert 2 galaxy were obtained due to the large 
field of view (FOV) of IBIS (19° x 19°, partially 
coded FOV) and SPI (35° x 35°, partially coded 
FOV), and because of the dithering observation 
strategy which is optimized for the spectrometer in 
order to allow a proper background determination. 
Observations were performed during seven INTE- 
GRAL revolutions in January, June and July 2003 
(Tab. 1). The total amount of exposure time is 512 
ksec with 274 ksec and 238 ksec in the January and 
June/ July campaign, respectively. Most of the ob- 
servations were carried out in a dithering mode. 
36 ksec of observation were carried out in staring 
mode, and are therefore only partly useful for SPI 
analysis, but they allow a more precise flux extrac- 
tion from the ISGRI data. Due to the smaller FOV 
of the two X-ray monitors ( JEM-X) and of the op- 
tical monitor (OMC) on-board INTEGRAL, NGC 
4388 has not been observed by those instruments. 
We will discuss the INTEGRAL data mainly in 
terms of the two observation campaigns in Jan- 
uary 2003 (JAN03) and Jimc/July 2003 (JUN03). 



2.1. IBIS 

The ISGRI spectrum has been extracted from 
ISGRI mosaics of the field created in 20 energy 
bands. The count rate and statistical error were 
extracted from the mosaic intensity and variance 
images for the ISGRI source position derived from 
the highest significant pixel of the mosaic in all en- 
ergy bands. The mosaic image fully redistributes 
the pixel values of the sky images of individual sci- 
ence windows that were created using version 3.0 
ot the Offline Science Analysis (OSA) software dis- 
tributed by the ISDC (Courvoisier et al. 2003b). 
Pixel values are weighted according to vigneting 
and variance. Finally we add a systematic error 
of 10% to the spectra. 

NGC 4388 was too faint to be seen by IBIS/PICsIT. 
Summing all available IBIS/ISGRI data together 
shows that NGC 4388 is detected up to 200 keV, 
though the data above 150 keV are consistent with 
zero flux on a la level. A flt to a single power law 
gives a photon index of T = 1.70^°°} (25 - 300 
keV) and a flux of 7.8 ± 0.5 mCrab in the 20 - 40 
keV energy band. Analysing the data according to 
the JAN03 and JUN03 campaign shows that the 
source was in a lower flux state in January 2003 
with a harder spectral slope (F = LBI^qJ^), while 
the June spectrum is steeper (F ~ 1.76 ^[JJ^) but 
has a higher flux. Note however that the spectral 
shape is still consistent with the F = 1.7 model on 
a Icr level. For details see Table 2. 

2.2. SPI 

In order to reduce the noise in the data reduc- 
tion process of the SPI data, only those pointings 
were considered where NGC 4388 was less than 15 
degrees off axis. This decreases the total amount 
of usable pointings to 122 with an exposure time 
of 178 ksec in JAN03 and 134 kscc in JUN03. OSA 
3.0 has been applied, except for some updates. In 
the binning of the data SPIHIST 3.1.2 was used 
in order to avoid a known problem with the ear- 
lier version. In the case of the image reconstruc- 
tion program SPIROS (Skinner & ConneU 2003), 
we used the more recent version 6.0.1, and also 
a response function which takes into account the 
results from the in-flight calibration on the Crab 
(instrument response function spi_irf_rsp_0015 and 
redistribution matrix spi_rmf_grp_0002). For de- 
tails of the SPI analysis procedures see Diehl et 



al. (2003). The SPI spectrum for the combined 
data of JAN03 and JUN03 is well described by a 
single power law with photon index F = 1.68 _„oc 
(20 - 500 keV) and a flux of 9 ± 2.6mCrab'in 
the 20 — 40 keV energy band, which is consistent 
with the ISGRI results. Analysing the JAN03 
and JUN03 subset it turned out that the source 
strength is not sufficient for SPI to achieve a rea- 
sonable signal to noise spectrum for the individ- 
ual epochs. In case of marginal detection the 
SPIROS software tends to overestimate the flux. 
This naturally leads also to a hardening of spectra, 
as the higher energy bins have a lower signal-to- 
noise. The results are shown in Table 2. It can be 
clearly seen that while the flux value of the com- 
bined spectrum is similar to that of ISGRI, the 
subsets with the lower signal-to-noise show higher 
flux values, and this trend is stronger at higher 
energies. Although the trend in flux is similar to 
the one observed by ISGRI, the spectral evolution 
is reversed. But the large error bars on the SPI 
results have to be taken into account and the re- 
sults are therefore still consistent with each other. 
For a detailed discussion of cross-calibration issues 
of the INTEGRAL instruments sec Lubihski et al. 
(2004). 

2.3. RXTE 

The All Sky Monitor (ASM) on board the Rossi 
X-ray Timing Explorer (RXTE) scans about 80% 
of the sky every orbit. This offers an unique way 
to monitor the emission of bright X-ray sources 
like NGC 4388 in the 1.5 ~ 12keV energy range. 
We extracted fluxes from the RXTE/ ASM data 
base. The fluxes have been averaged over the 
same time periods as the JAN03 and JUN03 IN- 
TEGRAL observations in order to have compara- 
ble results. The weighted mean of the JAN03 data 
was 3.5±1.2mCrab, while during the JUN03 cam- 
paign the flux was 2.2 ± 1.9mCrab. Due to low 
statistics the analysis in the three different ASM 
bands did not give further information about the 
flux variability. 

3. XMM-Newton observations 

Two observations with XMM-Newton (Jansen 
et al. 2001) were performed in July and Decem- 
ber 2002 as described in Table 3. Concerning the 
connection between the soft and hard X-ray emis- 



sion, the most interesting data for our study come 
from the two MOS cameras (Turner et al. 2001) 
and from the PN detector (Striider et al. 2001), 
as they cover the energy range 0.2 — lOkeV. The 
data have been reduced using the XMM-Newton 
Science Analysis Software version 5.4.1 and the 
model fitting using XSPEC 11.3 was done simul- 
taneously for the PN and MOS data. Both obser- 
vations are well represented by a Raymond-Smith 
model (Raymond & Smith 1977) at energies be- 
low ~ 2.5 keV with low abundance (7% solar) and 
a temperature of about 0.80 keV, applying only 
Galactic absorption (2.62 x 10^" cm^^). The so- 
called mekal model in XSPEC, developed by Mewe 
& Kaastra, gives the same results (T = 0.78 keV, 
7.4% Zq). 

Around 2.5 keV another additional component 
is seen in the data. This can be modeled by allow- 
ing the abundances of the single elements involved 
to vary. It turns out that this part of the spec- 
trum is best modeled when the sulfur abundance 
is about solar, thus being 20 times higher than the 
relative abundances of the other elements. 

The data do not allow to verify the more com- 
plex model of a photoionized gas applied by Iwa- 
sawa et al. (2003) to the Chandra data. This 
model avoids the (unrealistic) low abundance, and 
therefore seems to be more appropriate from a 
physical point of view. 

At harder X-rays, the spectrum is dominated by 
an absorbed power law. The column density of the 
cold absorber is 2.7 x lO^'^cm"^. The power-law 
index is not well constrained by the XMM-Newton 
data, so the photon index of P = 1.7 measured at 
hard X-rays, has been applied. In addition, a gaus- 
sian line at 6.39 ^p°J keV is apparent in both ob- 
servations. Details of the fit are listed in Table 4. 
This is consistent with a recent measurement by 
Chandra, which showed an Fe line of centroid in 
6.36 +°°2 keV (Iwasawa et al. 2003). Another 

weak line is detectable at 6.89^q'-^q kcV (probably 
Fe K/3) with a flux about ten times lower than the 
Fe Ka line. In both XMM spectra there was no 
sign of a further line at 7.1 keV, as reported from 
Chandra data by Iwasawa et al. (2003). However, 
due to the lower resolution of the XMM in this 
energy region, a contribution from a 7.1 keV line 
might be included in the wing of the Fe K/? line. 
The variation between the July and December 



data occurs mainly in the high energy part. This 
is clearly seen when comparing the December data 
with the model for the July data (Fig. 1). While 
there are little or no changes below ^ 2.5 keV, the 
spectrum varied significantly at higher energies. 
The iron fluorescence line at 6.39 keV was not af- 
fected by the variations. 

Sazonov & Revnivtsev (2004) studied the 
RXTE/FCA slew data of AGN in the 3-20 kcV 
energy band. This rcsuhed for NGG 4388 in 
a luminosity of L3_20koV — 4.5 x lO^^ergs"^, 
which is inside the range of that observed by 
XMM-Newton in 2002 (extrapolated L3-20koV = 
(2.8 - 7.0) X lO^'^ergs^^ for the July and Decem- 
ber 2002 observation, respectively). 

4. X-ray to Gamma-ray spectrum 

The INTECRAL data alone do not allow us to 
reconstruct the complete X-ray spectrum of NGC 
4388, mainly because the statistics of the ISGRI 
and SPI data do not constrain existing models. 

Comparison with former observations by SICMA 
(Lebrun et al. 1992), OSSE (Johnson et al. 1994), 
and BeppoSAX/VT)S (Risaliti 2002) shows a sim- 
ilar spectrum at hardest X-rays as observed with 
INTEGRAL. Measurements by BATSE gave a 
fluxof /20-iookoV = (2.6±1.5)xl0^^phcm^^s^^, 
compared with an INTECRAL value of (2.7 ± 
0.3) X 10~^phcm^^ s^^ This indicates that the 
soft gamma-ray spectrum does not vary dramat- 
ically as shown in the lightcurves for the 60 keV 
and 100 keV flux in Fig. 2. The large error in the 
SPI flux at 60 keV is mainly based on the fact 
that a strong detector background line complex is 
apparent in the 53 - 67 keV energy range resulting 
from ^"^Ge (Weidenspointner et al. 2003), which 
makes source flux extraction difflcult. The flux 
varied by a factor < 3 and the spectral shape is 
conserved during the 13 years of various observa- 
tions. 

In order to extend the spectrum into the en- 
ergy range below 20 keV we added the two XMM- 
Newton EPIC data sets, described in the previous 
section. 

The data were fit simultaneously in XSPEC 
11.3 applying the same model as for the XMM 
data alone: a Raymond-Smith model for a hot 
plasma with Galactic absorption, an absorbed 
power law, plus two Gaussian emission lines, rep- 



resenting the Ka and K/3 line. Figure 3 shows 
the combined spectrum of NGC 4388. The fit, 
which allowed for flux variations between the vari- 
ous non-simultaneous observations (see Sect. 5), 
gave a xt = 1-41 for f43f degrees of freedom. 
The normalization factors are in the range 2.5 
to 4.2 relative to the XMM-Newton July 2002 
data, which apparently represents the lowest flux 
level reported so far. The power law, dominat- 
ing the emission above ^ 2.5 keV, has a pho- 



ton index of F = 1.65 
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Using a more complicate model with a cut-off 
power law plus reflection from cold material (the 
so-called PEXRAV model; Magdziarz & Zdziarski 
1995) instead of the simple power law, does not 
improve the fit significantly. The data do not con- 
strain the power law and the reflection component: 
Fixing the photon index to F = 1.7 results in a 
folding energy of 380 ^^^ keV, indicating that a 
cut-off is possible but cannot be confirmed. 

Applying a single-power law plus an unresolved 
set of lines, represented by a broad Gaussian line 
around 1 keV, instead of the Raymond-Smith com- 
ponent does not give acceptable fit results {xl ^ 
2.0). 

Using only the high energy (> 20keV) data 
gives a single power law with F = 1.72||^gpg with 
xt — 1-5 (Fig. 4). Also in this energy region alone 
the more complex PEXRAV model does not give 
a significantly better fit. Statistics do not allow us 
to distinguish between a single power law and a 
PEXRAV model. The differences are smaller than 
the la error bars of the spectra as seen in Fig. 4 
and Fig. 3. 

5. Discussion 

The combination of the INTEGRAL data with 
previous observations by XMM-Newton, Bep- 
poSAX, CGRO, and SIGMA shows that the 
spectrum from 0.2 keV up to several hundred 
keV is well represented by a thermal Raymond- 
Smith model with only Galactic absorption ap- 
plied, plus a highly absorbed {Nh — 2.7 x 
10^^ cm^^) power law like emission (F = 1.7), 
and two Gaussian lines to model the iron flu- 
orescence line at 6.39 keV and 6.9 keV, respec- 
tively. The results are consistent with previous 
studies of NGC 4388 using subsets of these data 



(Risaliti 2002; Lebrun et al. 1992). Also obser- 
vations with Chandra (Iwasawa et al. 2003) and 
ASCA (Iwasawa et al. 1997; Forster et al. 1999) 
show a similar spectral behavior. 

The fact that the data of the different missions 
give an acceptable fit using the same normalization 
is based on the fact that the error bars in the high 
energy domain are still rather large. The accept- 
able correspondence between the ISGRI and SPI 
data leads to the assumption that the INTEGRAL 
calibration is reasonable and that normalization 
problems as reported e.g. in Courvoisier et al. 
(2003a) are not as important anymore. But large 
error bars on the high energy flux measurements 
(Fig. 2) still leave open the possibility of some er- 
rors in the calibration (Lubihski et al. 2004). 

The spectral data of NGC 4388 do not drop 
drastically above 200 keV, but are not of high 
enough significance to distinguish whether there is 
a cut-off, as expected from Comptonisation mod- 
els (Petrucci et al. 2001). On the other hand the 
absence of a cut-off would be in agreement with 
Deluit & Courvoisier (2003), who find a cut-off in 
the BeppoSAX/PDS spectra of Seyfert 1, but not 
in Seyfert 2. The XMM-Newton data show vari- 
ability only in the energy region above ^ 2.5 keV. 
This variability cannot be explained by a change 
in the absorption column density. But using the 
same model parameters of the July observation 
(Tab. 4) for the December 2002 data and allowing 
the normalization of the single power- law to vary, 
shows that a flux increase by a factor of 3.4 is suf- 
ficient to model the December data {xf, — 1.1; 
Fig. 1). The spectral shape of the power law, 
which describes the hard X-rays, the Fe Ka line 
strength, and the soft X-ray component do not 
seem to vary. 

The observed spectrum can be compared to 
the one derived by Ghisellini, Haardt, & Matt 
(1994) based on Monte-Carlo simulations for a 
Seyfert galaxy with inclination i = 60° — 63°, 
column density of 10^** cm^^ and a torus with 
half-opening angle of O = 30°. Based on these 
assumptions, an equivalent width of the iron fluo- 
rescence line of the order of EW ~ 100 eV would 
be expected, which would also mean a parallel 
evolution of the line flux with respect to the con- 
tinuum. But the Fe Ka does not show variability, 
while the underlying continuum varies by a factor 
of ^ 4 in the line region (see Fig. 5). The line 



flux of (6.7 ± 1.2) X 10-5phcm-2s-i and (7.8 ± 
1.4) X 10~^ phcm~^ s~^ for the July and December 
measurement is also consistent with the one mea- 
sured by Chandra ((9.3±1.9) x IQ-^phcm^^ g-i) 
(Iwasawa et al. 2003). The equivalent width 
therefore varied strongly {EW — 190 — 720 eV). 
The exception from this rule is the 1993 ASCA 
observation, which revealed a high continuum flux 
with simultaneous high iron line flux. But the 
rather stable line flux over several years with si- 
multaneous strong variations of the underlying 
continuum implies that the line emitting region is 
separated from the continuum source of the object 
by several light years. A similar behavior has been 
observed in i?X 77? spectra of some Seyfert 1 galax- 
ies, where the Fe Ka line also varies less strongly 
than the broadband continuum (Markowitz, Edel- 
son, & Vaughan 2003). 

Another difference of the spectrum observed 
here to the simulated one by Ghisellini, Haardt, 
& Matt (1994) is the fact that the spectrum of 
NGC 4388 does not seem to have a strong cut- 
off in the hard X-rays. This could be caused by 
a higher temperature of the corona (^ 100 keV) 
in the case of NGC 4388. If the non-detection of 
a reflection hump in the hard X-rays is real, this 
might be a hint for a non isotropic radiation or, 
generally speaking, a more complex geometry (see 
e.g. George & Fabian 1991). Apparently there is 
no scattering of the hard X-rays away from the line 
of sight, except for the iron Ka flourescence line. 
This indicates that in the case of NGC 4388 we in- 
deed observe the unscattered emission of the cen- 
tral engine at hard X-rays. Also no significant de- 
cline due to Klein-Nishina processes appears above 
50 keV can be detected in the present data. 

Most Seyfert galaxies show a softening of the 
X-ray continuum as sources brighten (Markowitz, 
Edelson, & Vaughan 2003; Nandra et al. 1997). 
The fact that the change in luminosity of NGC 
4388 in the hard X-rays, and therefore most prob- 
ably a change in the accretion rate of the central 
engine, is not accompanied by variations of the 
spectral shape remains enigmatic. One explana- 
tion might be a generally different radiation pro- 
cessing in Seyfert 1 and Seyfert 2, as indicated 
in e.g. Deluit & Courvoisier (2003). Matt et al. 
(2003) argue that absorption in Seyfert 2 galax- 
ies can be caused by several circumnuclear regions 
and favour a model with a possible temporary 



switching-off of the nuclear radiation to explain 
observed variability. Their model predicts an iron 
Ka equivalent width of ~ 100 eV for a Compton 
thick case with Nh — lO^'^cm"^, which is signif- 
cantly lower than the EW measured in NGC 4388 
with a similar absorption column density. 

The constant soft X-ray emission is most likely 
not linked to the AGN. This is also supported 
by the fact, that this emission is extended, as 
has been reported from ROSAT/URl measure- 
ments (Matt et al. 1994) and lately from Chan- 
dra (Iwasawa et al. 2003). The low abundance in 
the extended emission has been seen already in 
the ASCA observations, though with much larger 
uncertainties {Z = 0.05^q^^Zq; Iwasawa et al. 
1997). Iwasawa et al. (2003) argue that the ex- 
tended emission at low energies is more likely to 
origin from photoionized plasma. Even though it 
is not possible, based on the XMM-Newton data 
presented here alone, to give a preference to a ther- 
mal or photoionized plasma, the latter one has 
the advantage of avoiding the low abundance in 
the vicinity (within several kpc) of the AGN core. 
We therefore refer to the more simple model, ap- 
plied to the data presented here, which gives suffi- 
cient good fit results, even though it might be less 
meaningful in a physical context. The Chandra 
data also allowed a space resolved study at soft 
X-rays, which was not possible with XMM^ and 
found that the AGN core itself apparently con- 
tributes less than 10% to the soft X-ray flux (0.5 
- 2.0 keV) compared to the surrounding gas. 

The INTEGRAL /I'&GRl measurements indi- 
cate that a lower flux in the soft gamma-ray region 
is accompanied by a harder spectrum, but the er- 
rors on the spectral slope measurement are too 
large to be statistically significant. Since the SPI 
data at the same time show the contrary spectral 
evolution, it is likely that the spectral slope did 
not change and that there is only a change in flux 
(20 - 200 keV) by a factor of - 1.5. 

This variation is not seen in the RXTE/kS>yi 
data, which cover the softer energy range (1.5 - 
12 keV). Note also the large error values, which 
are consistent with no flux variation for both, the 
INTEGRAL and the RXTE data. 



6. Conclusion 

The complex X-ray spectrum of NGC 4388 is 
composed out of three major components, which 
originate from different regions: 

• Hot plasma component: Emission be- 
low ^ 2.5 keV is dominated by emission of 
a hot (10^ °K), optically thin plasma with 
low abundance (about 7 % solar) . This com- 
ponent does not seem to be variable. In 
ROSAT/RRl and Chandra observations this 
emission appeared to be an extended X-ray 
nebula out to several kpc from the AGN. 
Chandra observations imply that this com- 
ponent might be due to photoionized gas. 

• Hard X-ray emission: Emission above 
~ 2.5 keV follows a highly absorbed (Nh ~ 
2.7 X 10^^ cm^^) power law with photon in- 
dex r = 1.7 and is related directly to the 
AGN. While the flux of this component var- 
ied by a factor ~ 4 over the past ten years, 
the spectral shape has been constant within 
the error bars. No reflection component or 
cut-off is detectable in the data so far. 

• Iron line emission: The Fe Ka fluores- 
cence line at 6.39 keV and K/3 at 6.9 keV 
are both consistent with the redshift of the 
AGN. The former had a constant flux of 
~ 7.5 X lO-'^ ph cm-2 s^^ from 1995 to 2002. 
The iron line emission is decoupled from the 
broadband continuum, similar to the behav- 
ior observed in several Seyfert 1 galaxies 
(Markowitz, Edelson, & Vaughan 2003) and 
NELG(e.g. NGC 5506; Bianchi et al. 2003). 
Also the fact that no line broadening is seen 
supports that it is probably generated sev- 
eral light years away from the AGN core, for 
example in the molecular torus. 

The INTEGRAL data do not allow to constrain 
the model on the high energy emission, i.e. there 
is no sign (except the iron fluorescence line) of 
a Compton reflection component or a cutoff at 
higher energies but its existence can also not be 
ruled out. 

Following the observations reported here, NGC 
4388 was studied in a dedicated INTEGRAL ob- 
servation, which was performed in staring mode in 



July 2003. Though staring mode observations usu- 
ally do not produce useful SPI spectra for sources 
as faint as NGC 4388, the combination of ISGRI 
and JEM-X data might show in more detail the 
connection between the hard and soft X-rays, as 
seen in this work from simultaneous INTEGRAL 
and RXTE data. Another 100 ksec TOO observa- 
tion in June 2004 and 500 ksec open time obser- 
vation of 3C273 in INTEGRAL'S AO-2 might give 
further insights into the cut-off of the NGC 4388 
spectrum at high energies, e.g. when combined 
with the already existing data presented here. 

This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory. We thank 
the anonymous referee for constructive comments, 
which helped to improve the manuscript. 
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Table 1: INTEGRAL observations 



obs. date 


INTEGRAL 


SPI 


number of 


observation 


start 


revolution 


exp. time [sec] 


pointings 


mode 


05/01/03 


28 


80520 




39 


5x5 


11/01/03 


30 


10558 




2 


staring 


17/01/03 


32 


182643 




43 


5x5 


16/06/03 


82 


20000 




1 


staring 


06/07/03 


89 


152950 




47 


5x5 


09/07/03 


90 


59989 




17 


5x5 


18/07/03 


93 


5283 




2 


staring 



Table 2: Fit results for a single-power law model 



Observation 


Photon 


F% 


F^x 


Fl 


Xiidof) 




Index 


20-40 keV 


40-100 kcV 


100-200 keV 




Jan.03 / ISGRI 


1 61 +'^-^'^ 

-^•"^ -0.18 


4.5 ±0.9 


8.2 ±3.1 


8.5 ±5.9 


1.16 (8) 


Jan.03 / SPI 


1 56 +°'^^ 
-■^•^0 _Q53 


9.7 ±3.6 


17.0 ±7.8 


19.9 ±6.2 


0.78 (2) 


Jun.03 / ISGRI 


1 76 +°-^^ 


6.6 ±1.2 


10.6 ±2.8 


9.7 ±4.5 


1.42 (8) 


Jun.03 / SPI 


1 00 +0-37 
^■^^ -0.27 


11.0±3.0 


24.3 ±8.0 


34.0 ±11.6 


0.70 (5) 


summed / ISGRI 


1.70 tTo\ 


7.0 ±0.4 


11.7±1.4 


11.3 ±3.3 


1.69 (14) 


summed /SPI 


1-68 t',% 


7.7 ±2.1 


12.2 ±4.8 


12.6 ±7.7 


1.30 (4) 



"■ un-absorbed flux in 10 ^^ erg cm 



Table 3: XMM-Newton observations 



obs. date 
start 



exp. 



MOS-1 
time [seel 



MOS-2 
exp. time [sec] 



PN 

exp. time [sec] 



07/07/02 
12/12/02 



9850 
11667 



9937 
11667 



4540 
8291 



Table 4: Fit results for the XMM-Newton EPIC data 
Observation T[kT] abundance Nh 



+0.04 

-0.05 
+0.03 

-0.03 



+0.03 
-0.02 
+0.01 
-0.01 



[10 



23, 



+0.20 
-0.21 
+0.07 
-0.07 



[10 



/fcKq 

'"^ phcm" 



EW 

[eV] 



a 

[eV] 



Xlidof) 



July 2002 
Dec. 2002 



0.81 
0.80 



0.08 
0.06 



2.45 
2.79 



6.7 ±1.2 566 67 ±16 1.06(526) 

7.8 ±1.4 220 69 ±12 1.25(861) 



NGC 4388 combined spectrum 
XMM-Newton EPIC MOS/PN, BeppoSAX/PDS. INTEGRAL SPI/ISGRI, OSSE, and SIGMA 



NGC 4388 XMM-Newton spectrL 
December 2002 dota versus July 2002 model 




channel energy (keV) 



Fig. 1.— Comparison of XMM-Newton MOS data 
from December 2002 with respect to the model for 
the July 2002 EPIC/MOS data. 




1990 



1995 2000 

time [year] 



Fig. 2.— Lightcurve for NGC 4388 at 60 keV (up- 
per solid line) and at 100 keV (lower dashed line). 
Data taken from SIGMA (circle), OSSE (pen- 
tagon), BeppoSAXlVDS (hexagon), SPI (square), 
and ISGRI (triangle). 




Fig. 3.— NGC 4388 combined spectrum {E'^Je 
vs. E), including XMM-Newton Y^VIC MOS/PN, 
BeppoSAX PDS, CGRO OSSE, SIGMA, and IN- 
TEGRAL ISGRI/SPI. The high energies are dom- 
inated by the highly absorbed reflection compo- 
nent plus the Ka and K/3 iron fluorescence lines, 
while the lower energies show emission character- 
istic for a hot plasma with low abundance. XMM 
data from July 2002 are not shown in this plot for 
better visibility. 



NGC 4388 high en&rgy spectrL 
BeppoSAX/PDS, INTEGRAL ISGRI/SPh OSSE, ond SIGMA 




100 200 

annel energy (keV) 



Fig. 4. — NGC 4388 combined high energy 
spectrum, including BeppoSAX PDS, GGRO 
OSSE, SIGMA, and INTEGRAL ISGRI/SPI. 
The spectrum is shown in instrument dependent 
counts s^^ keV^"'^. 
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Fig. 5. — Flux in the 2-10 keV band (triangles) 
and of the 6.4 kcV Fe Ka line (circles) measured 
by ASCA (1993 and 1995), BeppoSAX (1999 and 
2000), Chandra (2001), and XMM-Newton (2002). 
The underlying continuum is highly variable, while 
the line flux shows no significant variation, except 
for the ASCA measurement in 1993. 
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